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Abstract 

We consider two dilute gas Bose-Einstein condensates with opposite velocities 
from which a monochromatic light field detuned far from the resonance of the 
optical transition is coherently scattered. In the thermodynamic limit, when 
the relative fluctuations of the atom number difference between the two con- 
densates vanish, the relative phase between the Bose-Einstein condensates 
may be established in a superposition state by detections of spontaneously 
scattered photons, even though the condensates have initially well-defined 
atom numbers. For a finite system, stochastic simulations show that the mea- 
surements of the scattered photons lead to a randomly drifting relative phase 
and drive the condensates into entangled superpositions of number states. 
This is because according to Bose-Einstein statistics the scattering to an al- 
ready occupied state is enhanced. 
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I. INTRODUCTION 



Bose-Einstein condensates (BECs) of ultra-cold trapped atomic gases [0-^] have stim- 
ulated interest on the coherence properties of matter. BECs are expected to exhibit a 
macroscopic quantum coherence which in thermal atomic ensembles is absent. Even BECs 
with a well-defined number of atoms, and with no phase information, could show phase 
correlations in particular measurement processes on atoms P-[TT[], or on photons The 
relative phase between two BECs could be determined, for instance, by various optical meth- 
ods [P^-fTB|]. In recent experiments Andrews et al. JDJ have found evidence of macroscopic 
quantum coherence in a BEC by measurements of the interference of two condensates by 
absorption imaging. The two independent and spatially separated BECs were created by a 
repulsive optical force in the center of the trap. 

In this paper we consider an optical analogue of the atom detection schemes of two BECs 
in different momentum states • Instead of looking at the spatial interference pattern 

we combine the scattered photons from the atomic transitions between different BECs with 
a photon beamsplitter. The measurements on scattered light has evident advantages over 
atom counting from a theoretical point of view. In the case of light scattering we can use 
the well-known theories of photon detection ||18| . Also, the measurement of spontaneously 
scattered photons is nondestructive for the condensates, because only light is scattered and 
no atoms are removed from the two BECs. In our measurement scheme it is shown via 
the simulations of stochastic Schrodinger equations that the detections of spontaneously 
scattered photons drive the condensates into macroscopic quantum superpositions of phase 
and number states ("Schrodinger cats"). The phase superpositions are a consequence of the 
particular measurement process, which is insensitive to certain phase values. The entangled 
number state superpositions follow from the properties of Bose-Einstein statistics and from 
the macroscopic quantum coherence of BECs. The number state superpositions are multi- 
particle quantum states with spatially nonlocal correlations. Due to the large fluctuations 
of the number difference between the BECs the relative phase drifts randomly in the case 
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of a finite system and no stable phase is built up by measurements. The fluctuations of 
the number difference vanish in the thermodynamic limit and the measurements establish a 
stable phase. 

Recently, Cirac et al. have studied the ground state of two coupled BECs by varia- 



tional techniques. They have found that under certain conditions the state with a minimum 
energy corresponds to a macroscopic superposition of number states. 

We begin in Sec. [II A| by introducing basic relations. In the limit of large detuning 



of the driving light field from the atomic resonance the excited state operators may be 
eliminated adiabatically. We obtain an effective two-state Hamiltonian coupling the two 
BECs. We study the dynamics of the system in terms of stochastic trajectories of state 
vectors. In Sec. [II B| we consider the thermodynamic limit, where the fluctuations of the 
number difference between the BECs vanish. The results of simulations for a finite system 
are presented in Sec. |II C[ In Sec. [II D| we show that the number state superpositions could 



be detected by considering the intensity correlations of the scattered light. Finally, a few 
concluding remarks are made in Sec. fTJ. 

II. OPEN SYSTEM DYNAMICS 
A. Basic relations 

The internal quantum state for both condensates is denoted by \g). This state is optically 
coupled to the electronically excited state \e) by the driving electric field with a dominant 
frequency Q. The light field is assumed to be in a coherent state and detuned far from the 



resonance of the atomic transition. The two BECs are assumed to be optically thin p0[ and 
in the momentum states ko and — ko. We consider the situation in which the condensates 
are overlapping when the light is switched on. We only consider the coherent spontaneous 
scattering between the condensates, which is stimulated by a large number of atoms in the 
condensates. By spontaneous scattering we mean that the emission is not stimulated by light, 
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although it is stimulated by atoms. The decay into non-condensate center-of-mass (cm.) 
states is also stimulated by the Bose-Einstein statistics. However, at very low temperatures 
this stimulation is much weaker because most of the particles are in the condensates. In 
addition to the Bose stimulation of spontaneous emission there is unstimulated free-space 
decay, at rate 7, which is always present. With a sufficiently large number of atoms in the 
two BECs the free-space decay may be ignored. 

The annihilation operators for the two BECs are g> ko and <7_ ko . Here g ko denotes the 
annihilation operator for the electronic ground state \g) and the cm. state ko with the 
corresponding wave function ff)ko (r). To simplify the notation we define b = # ko , 0&(r) = 
0g,k o ( r ), c = g_ ko , and c (r) = <p g _ ko (r). We obtain for the Hamiltonian [pC|[2T 

H = hei Q tfb + he 9 _ ko c ] c + £ h{u eg + e k ) e k e k + ]T %u q a\ 

k 

- £ ( />r d 9e ■ E(r)0*(r)0 ek (r) tfe k + H.c. 

k \ J 

- E (/ d3r d 9e ■ E(r)^(r)0 ek (r) Je h + H.c.) , (1) 

where the excited state wave function for the cm. state k is ek . The dispersion relations for 
the ground state and excited state cm. frequencies are e k and e k , respectively. The photon 
annihilation operator for the mode q is a q . The internal atomic energy is described by the 
frequency u eg of the optical transition between the electronic ground state and excited state. 
The last two terms in Eq. (0) are for the atom-light dipole interaction. The dipole matrix 
element for the atomic transition e — >• g is given by d ge . We consider the translationally 
invariant system, where the eigenfunctions for the condensates are plane waves: </>&(r) = 
e lk °' r /\/V and c (r) = e~ lk °' r / y/V . The driving electric field is also described by a plane 
wave Ej(r) = Se e' ; ( k - r -™)/2. 

In the limit of large detuning, A = Q — u eg , the excited state operators e k in Eq. (P 
may be eliminated adiabatically, and the cm. energies of the excited state may be ignored 
The system may then be described by an effective two-state Hamiltonian: 

H = hel o tfb + he 9 _ ko c^c + ^ a \ a q 



d " rdae ' E(r) de9 ' E ( r )^(r)<Mr) 
+ (tfc J d 3 rd 9e -E(r)d e9 -E(r)0*(r)0 c (r) + H.c.)| . (2) 

Here we have used the fact that for the plane waves 0£(r)0b(r) = 0*(r)0 c (r). The total 
atom number operator is given by N = tfb + c^c. Because the total number is conserved, 
the operator N contributes to the measurements only through a constant phase shift. Thus, 
we may ignore the term proportional to N in Eq. (Q). 

We consider the dynamics of the two BECs and the driving light field as an open quan- 
tum system and eliminate the vacuum electromagnetic fields. The set-up of our Gedanken 
experiment is given in Fig. |l|. The incoming light field is scattered from two overlapping 
BECs moving with opposite velocities. The scattering processes in which an atom scatters 
back to the same condensate introduce the term proportional to the total number of atoms 
in Eq. and they may be ignored. In the scattering processes in which atoms scatter 
between different condensates the light beams are deflected due to the recoil momentum. 
In Fig. [I] a photon deflected to left corresponds to the change of the momentum of an atom 
from — ko to ko, i.e. the amplitude of the scattered electric field is proportional to b*c. 
Similarly, a photon deflected to right corresponds to the change of the momentum of an 
atom upon scattering from k to — k . In this case the amplitude of the scattered electric 
field is proportional to c'b. The scattered light beams are combined by perfectly reflecting 
mirrors and a 50-50 photon beamsplitter. The detection rate of photons on the detectors is 
the intensity of the scattered light J(r) = 2ceo(E _ (r) • E + (r)) integrated over the scatter- 
ing directions divided by the energy of a photon tick. Writing the electric fields in the far 
radiation zone (kr ^> 1) |2(J we obtain the detection rate at the channel j: 



lj ~~ Tick 



l -j ~dQ il r 2 I j (r) = 2T{C}C j ), (3) 



r = ^ id g\ 2 (4) 



where the linewidth of the electric dipole transition is given by 7 = d 2 g k 3 / (6irheo). The two 
relaxation channels corresponding to the two output channels of the beamsplitter are 



C\ — J x , C*2 — Jy , (5) 

where the familiar angular momentum operators obeying SU(2) algebra are defined by 

X = l(tfc + Jb) , (6a) 

J, = 1 (V c - c tft) , (6b) 

■/^(fetfe-c+c) , (6c) 

and J 2 = J 2 + J 2 + J 2 = (N/2 + l)N/2 is the Casimir invariant. 

The system Hamiltonian for the BECs and for the driving electromagnetic field is elim- 



inated completely in the interaction representation and according to Ref. |12 we may then 
write down the equation of motion for the reduced density matrix of the system in the limit 
of large detuning of the driving light from the atomic resonance 

Ps = -r E (cjQps + psCjQ - 2C iPs Cf) . (7) 

i=\ 

If we assume the condensates to be in coherent states with equal mean atom numbers 



(b) = ^N/2e tipb and (c) = yN/2e l<Pc , the intensities of the scattered light in the two 
channels are 

Ji oc (J 2 ) oc (cosy?) 2 , J 2 oc (J 2 ) oc (simp) 2 . (8) 

Here we have defined the value of the relative phase by (p = ip c — ip^, where ipi and (p c are 
the macroscopic phases of the condensates b and c, respectively. There is an ambiguity in 
Eq. (H) between the phase values ±ip and tt ± ip. For phase-sensitive homodyne detection 
this ambiguity vanishes. 

The dynamics of the density operator from Eq. ([?[) may be unraveled into stochastic tra- 
jectories of state vectors p2|-[2"4|j. The procedure consists of the evolution of the system with 
a non-Hermitian Hamiltonian H c g, and randomly decided quantum 'jumps' corresponding 
to the direct detections of spontaneously emitted photons. The system evolution is thus 
conditioned on the outcome of a measurement. The non-Hermitian Hamiltonian is obtained 
from Eq. (0) 
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ihV J2 C}Cj = -i%Y ( J 2 - J 2 2 ) 



(9) 



The non-Hermit ian Hamiltonian H e g determines the evolution of the state vector ip sys (t). 
If the wave function ip S y S (t) is normalized, the probability that a photon from the output 
channel j (j = 1, 2) of the beamsplitter is detected during the time interval [t, t + 5t] is 



The probability of no detections is 1 — P\ — Pi- 

The implementation of the simulation algorithm is similar to Ref. [|l^]. At the time 
to we generate a quasi-random number e which is uniformly distributed between and 1. 
We assume that the state vector ip sys (to) at the time to is normalized. Then we evolve 
the state vector by the non-Hermitian Hamiltonian H e R iteratively for finite time steps 
At ~ 5t. At each time step n we compare e to the reduced norm of the wave function, 
until (ip S y S (t + nAt) \ ip sys {t + nAt)) < e, when the detection of a photon occurs. After the 
detection we generate a new quasi-random number rj. We evaluate Pi and P2 from Eq. ( |T0D 
at the time of the detection. If 77 < P\j [P\ + P 2 ) we say the photon has been detected from 
channel 1. If the photon has been observed during the time step t — > t + At we take the 
new wave function at t + At to be 



which is then normalized. Otherwise, 77 > P\j \P\ + P2) and the photon has been detected 
from channel 2. In that case the new wave function before the normalization reads 



Pj{t) = 2r(Ws(t) I qcj 1 </> sys (t)) st . 



(10) 




(ii) 




(12) 



After each detection the process starts again from the beginning. 



B. Thermodynamic limit 



Before presenting numerical results of the simulations of the stochastic Schrodinger equa- 
tions, we investigate qualitatively the build-up of the macroscopic coherence by the mea- 
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surement process in the limit iV — > oo. The eigenstates of the effective Hamiltonian H e g in 
Eq. (^|) are number states which have flat phase amplitudes. Thus, the time evolution of 
H c ff does not support any particular phase value over other values, and the relative phase 
between the two BECs is determined by the distribution of the photon detections between 
the two output channels of the beamsplitter. Because the two relaxation channels do not 
commute, [J x , J y ] = i.J z , the state of the system depends also on the particular order in 
which the scattered photons are detected. This complicates the analysis substantially. An 
evident simplification is to consider the thermodynamic limit N — > oo, where the relative 
fluctuations of the number difference between the BECs vanish (J z ) / (N) — ► 0. In the next 
section we consider in the numerical simulations systems which are far away from this limit. 
However, it turns out that the qualitative behaviour of the macroscopic phases is still very 
similar. 

In the thermodynamic limit we can replace the angular momentum operators J x and J y 
by (N/2) cosy? and (N/2) simp respectively, where is the relative phase operator between 
the two BECs. Then the two relaxation channels in Eq. (|5|) commute, [cos if, sin if] = 
AiAN /N 2 — > 0, and the relevant commutation relations are given by 

[AN, cos <f] = i sin <p, [AN, sin ip] = — i cos (p , (13) 

where we have written AN = J z . Geometrically, the thermodynamic limit may be under- 
stood as a restriction of the dynamics of the angular momentum variables to the equator 
(J z ) ~ of the Bloch sphere, as the radius of the sphere goes to infinity. 

Now we can use the procedure developed in Ref . [[Uj . To simplify the notation, we ignore 
the spatial dependence of the wave functions. We expand the number state \N/2,N/2) in 
terms of the overcomplete set of phase states ||25|| : 



W)n = (& f e-^ 2 + cte^T |0> . (14) 



A similar analysis to Ref. [TO] then leads to the state of the system after n\ and ri2 detections 



from output channels 1 and 2 of the beamsplitter respectively: 
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t-2-K 

\i){n u n 2 )) oc / (cosy?)™ 1 (sin <^) n2 \(p) N . (15) 
Jo 

The value of the phase ip that maximizes the integrand satisfies the relation tan 2 ip = n 2 /ni. 
If < ipo < 7r/2 is a solution for the maximum amplitude, then — ip and 7r ± y? are also 
solutions. For ni,ri2 3> 1, we can express the integrand in Eq. ( |T5| ) in terms of exponential 
functions and expand the exponents in a Taylor series around ±<p and n ± ip . We obtain 

|V(ril,n 2 )) OC ^ dip { e -^-W>) 2 + (.^n^-n^-w) 2 

+ (_l)«i+™2 e -n(v+— ^o) 2 + (_ 1 )" 2e -n(^ ) 2 | ^ ^ ( 16 ) 

where n = rii + n 2 . The phase distributions are Gaussians centered at four superposition 
values ±<y9 and tt ± (p , where < (po < tt/2 is a solution for tan 2 ip = n 2 /n 1 . The phase 
is well-defined with a narrow width for the Gaussian. The transition from the binomial 
distribution of Eq. (|15|) to the normal distribution of Eq. ( |16|) is just the realization of the 
Central Limit Theorem for a large number of detections, while the superpositions are a 
consequence of the particular detection method, which is insensitive to the phase values ±ip 
and Ti ± ip according to Eq. (^) . 

C. Numerical results 

For a finite system the two relaxation channels from Eq. (||) do not commute and the state 
of the BECs depends on the particular order in which photons from the two output channels 
of the beamsplitter are detected. We have simulated the measurements of the spontaneously 
scattered photons numerically for N = 200 atoms. Even though we start from the initial 
number state Nb = N c = 100 with no phase information, the detections establish coherence 
properties for BECs similarly to Eq. (|T6|). However, the value of the phase (po in Eq. (^) does 
not stabilize due to the moderate value of iV chosen, even for a large number of detections. 
The two BECs are also in entangled number state superpositions. The coherence properties 
of the BECs vary strongly even during single realizations of the measurement process. In 
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the extreme case the condensates approach an entangled number state with almost all the 
atoms in one of the two BECs. 

The emergence of the number state superpositions may be understood from the quan- 
tum statistical properties of Bose-Einstein particles. Because the scattering to the non- 
condensate modes is ignored, the total number of atoms in the two BECs is conserved and the 
atom numbers are entangled. According to Bose-Einstein statistics the scattering to an al- 
ready occupied state is enhanced. For the initial state | N/2, N/2 ) the probability for light to 
scatter atoms between the two BECs is, in the case of spatially overlapping condensates and 
in the limit of large number of atoms, approximately proportional to (N/2) 2 . Because the 
detected photons corresponding to the two atomic transitions between the BECs are indistin- 
guishable, the number state distribution remains symmetric with respect to the initial state 
during the scattering process. For an entangled number state (\ N — k, k) + \ k, N — k)) / \J2 
the scattering probability is approximately proportional to (N — k)k < (N/2) 2 . Hence the 
states with unequal atom numbers have smaller scattering rates and they are more stable. 
It should be pointed out that it is not necessary to have initially equal atom numbers to 
obtain number state superpositions, although the distributions are perfectly symmetric only 
if the initial atom numbers are the same. 

The state of the BECs may be described in terms of quasiprobability functions. For the 
number state distribution of atoms \ipb) = J2n c n\ n ) i n the condensate b we have evaluated 



the Q function [26 



Q(a) = 



J N a n c* 2 

. /TT 



(17) 

In Figs. H, |3], and f|, we have plotted |^&|, the absolute value of the wave function in the 
condensate b in the number state basis, and the corresponding Q function at different times 
during a single realization of measurements. Figures |2] and |3] represent typically observed 
distributions, when several thousands of detections are made. In Fig. [| we have a special 
case in which almost all the atoms are in one of the two BECs. In Fig. two distinct 
peaks in the number distribution are clearly observed. The first peak is centered at Nb ~ 30 
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atoms, i.e. N c ~ 170 atoms, and the second at Nb — 170 atoms, i.e. N c ~ 30 atoms. Only 
odd number states are occupied. This is because in each photon detection the states with 
only even numbers in the atom number distribution are changed to the states with only 
odd atom numbers and vice versa. In particular, for a coherent system we may define even 
| a, +) and odd | a, — ) coherent states by | a, ±) oc | a) ± | — a) ]27J. These are states which 
have only even or odd numbers in the atom number distribution and they correspond to 
superposition states with two different phase values shifted by ir. 

In Fig. §d we have plotted the corresponding Q function from Eq. (fL7|). The Q function 
gives the phase-space distribution. The amplitude and phase quadratures are denoted by X 

1/2 

and Y . In polar coordinates the radius in the xy plane is equal to N b and the polar angle 
is the relative phase between the two BECs. In Fig. it is easy to see the two different 
sets of peaks corresponding to the two dominating values in the number distribution. All 
the peaks are aligned parallel to the x axis. This is the reason that two of the four different 
phase values from Eq. ([16]) are indistinguishable. Although the distribution in Fig. ^|a is 
symmetric, the number squeezing of the peak with the larger atom number is much stronger. 
The fringes indicating a quantum interference in the Wigner distributions [2B| are absent in 



the Q representation, so that graphs of Q functions do not obviously distinguish between 
pure states and statistical mixtures. However, because we are dealing with basis vectors 
instead of with density matrices, it is evident that we have a pure state. 

In Fig. |^ the distribution in the number state basis and the corresponding Q function 
are plotted in the same run of measurements as in Fig. ^ but at different time. In Fig. ^a 
two distinct peaks in the number distribution are not as far apart as in Fig. 0a. In the 
Q representation, in Fig. [3]b, it is easy to see that the value of the relative phase between 
the BECs is different from Fig. [| In Fig. |3]b all the four phase values from Eq. are 
clearly observed. The value of the relative phase between the condensates wanders during 
the simulations and does not stabilize to any definite value. In Fig. [| we have one more 
graph from the same run of measurements. In this case the BECs are almost in an entangled 
number state with all the atoms in one of the two condensates. Because the state of the BECs 
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is closer to a number state than to a coherent state, the relative phase is not well-defined. 

In the calculations we only considered the atom-stimulated scattering to the BECs. It is 
not necessarily a well-justified assumption to ignore the unstimulated free-space decay if the 
condensates contain only 200 atoms. However, the purpose of the numerical simulations was 
to demonstrate the general properties of the finite systems with a convenient computational 
efficiency In the simulations the physical behaviour remained qualitatively the same even 
though the number of atoms was significantly increased. In the real experiments BECs have 
typically contained many more than 200 atoms 

The interaction of the BECs with their enviroment creates dissipation and the deco- 



herence of the macroscopic superpositions [pq| . Decoherence by amplitude damping or by 



phase damping has been estimated by Walls and Milburn p9| . The amplitude damping 
corresponds to the losses of atoms from the BECs. In this case the off-diagonal elements 
of the density matrix between two coherent states may be shown to be dephased by the 
factor (a|/?) 1 ~ cxp (~ A '), where A is the loss rate for atoms. The phase damping may, e.g., 
be a consequence of elastic two-body collisions. In this case the off-diagonal elements of 
the density matrix between two coherent states with unequal atom numbers Ni and N 2 are 
damped by the factor exp {— \{N\ — N 2 ) 2 t/2}. These are examples of the decoherence of the 
ensemble averages over the measurement processes. If decoherence can be associated with 
measurements, evolution of single realizations may be analyzed by stochastic evolutions of 



state vectors IpOf . Although the decoherence of number state superpositions of the BECs 



may exhibit some interesting features, we do not consider this in the present paper. 

D. Detection of number state superpositions 

In this section we consider the detection of the number state superpositions. The phase 
superpositions could in principle be measured by simply interfering the condensates. How- 
ever, the different phase superpositions correspond to either even or odd coherent states. 
As explained previously, these are states which have only even or odd numbers in the atom 
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number distribution. Thus, in practice the losses of atoms from the BECs could shift the 
fringes and wipe out the qualitative features from the interference pattern. 

The existence of the number state superpositions could be verified, for instance, by 
considering the intensity correlations of the scattered light from the two BECs. If we assume 
that the condensates are flying apart and that they are already spatially separated, the 
amplitude of the spontaneously scattered light field from the condensate b has roughly 
the dependence |Ejj~| oc £d eg /(hA) b*b, and from the condensate c |E+| oc £d eg /(hA) c^c 
20fl . Here we have again considered only the coherent spontaneous scattering of atoms to 



the BECs stimulated by large atom numbers, i.e. scattering to the non- condensate cm. 
states has been ignored. Because the BECs are now spatially separated, only the scattering 
processes in which an atom scatters back to the same condensate are included. 

For a number state | N/2, N/2 ) with large N, the intensities of the scattered light from 
the spatially separated condensates b and c are approximately (/&) ~ (I c ) oc (N/2) 2 . The 
intensity correlations satisfy (hie) (N/2) 4 . For an entangled number state (\N — k,k) + 
| k,N-k))/\/2v?e have (I b ) ~ (J c ) oc {(N- k) 2 + k 2 }/2 > (N/2) 2 and (I b I c ) oc (N-k) 2 k 2 < 
(N/2) 4 . An especially interesting case is the situation in which the superpositions are far 
apart: k < N/2. Then, {(N - k) 2 + k 2 }/2 ~ 2(N/2) 2 and (N - kfk 2 < (iV/2) 4 . Thus, 
for the present case of spatially separated condensates the scattering rate from each BEC 
is larger than the scattering rate given by N/2 atoms; on the other hand, the intensity 
correlations are at the same time strongly reduced. These conclusions are also valid for the 
case where the superpositions are of coherent states instead of number states, as long as the 
overlap between the superpositions is negligible. 

Because the detection of the number state superpositions relies on atom stimulated 
scattering to the BECs, the entanglement between the condensates is not destroyed in the 
measurement process. However, the light scattering still creates decoherence by phase damp- 
ing explained in the previous section. This decoherence may be reduced by balancing the 
detection rates of scattered light from the BECs. 
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III. CONCLUSIONS 



We have shown that two BECs can be driven into macroscopic superpositions of num- 
ber and phase states by measurements of spontaneously scattered light. The number state 
superpositions are entangled and spatially nonlocal "Schrodinger cat" states with high oc- 
cupation numbers. No stable relative phase between the BECs is established for a finite 
system and in the extreme case the condensates approach a number state with almost all 
the atoms in one of the BECs. This is an example of the strong effect of measurements on 
the state of the condensates. For a finite system detections necessarily perturb the phase 
and it is not irrelevant what particular phase measurement process is used; because different 
measurement procedures may affect the system in a very different way, it is not evident a 
priori what kind of coherence properties, if any, are established in a detection process. 



In the system considered in Ref. [12] two BECs are in two different Zeeman levels and two 
phase coherent laser beams drive Raman transitions between the condensates. In that case 
the relative phase between the two BECs is established by measurements of spontaneously 
scattered photons, even though the condensates have initially well-defined numbers of atoms. 
In the present paper the large fluctuations of the number difference between the BECs in the 
case of a finite system lead to a randomly drifting relative phase. The measurement scheme 
considered here, with only one laser beam and the two BECs differing in their external 
quantum numbers, is closer to the experimental set-up used at MIT in a nondestructive 
optical detection of a BEC [31 , 17| . In the phase-contrast [p/?]] or dark-ground [31] imaging 
of the BECs the role of the mirrors is played by a lens. Although so far all measurements 
of the interference pattern of BECs have been destructive, nondestructive measurements 
could possibly be performed in the near future. Only with nondestructive imaging could 
one measure how the system evolves in time as a result of the detection process. 
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FIGURES 

FIG. 1. The experimental set-up. The incoming light field is scattered from two overlapping 
BECs moving with opposite velocities. The atoms scattering from one condensate to another 
change the momenta of the scattered photons. The scattered photons are collected by reflective 
mirrors and a 50-50 beamsplitter. The photons are detected from the two output channels of the 
beamsplitter. The photons scattered forward introduce only a constant phase shift and they are 
ignored. 

FIG. 2. Stochastic simulations of the detections of spontaneously scattered photons for 200 
atoms. A typical distribution of (a) the absolute value of the wave function |^| in the condensate 
b in the number state basis, and (b) the corresponding Q function after approximately 5000 detec- 
tions. In (a) two distinct peaks in the number distribution correspond to entangled number state 
superpositions. The peaks are centered at A5 ~ 30 and at iV& ~ 170 atoms. In the phase-space 
plotting of the Q function (b) the radius in the xy plane in the polar coordinates is equal to N^ 2 
and the polar angle is the relative phase between the two BECs. The four peaks correspond to the 
two dominant occupation numbers and two different phase values. 

FIG. 3. Another representative graph from the same run of measurements with (a) the number 
state distribution and (b) the Q function. The value of the relative phase has changed from the 
previous figure. The two entangled number state superpositions and all the four phase values are 
clearly observed in the Q representation. 

FIG. 4. The same run of measurements as in previous figures. The plotting of an extreme case 
in which the BECs are in an entangled number state with almost all the atoms in one of the two 
condensates. 
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Fig.l Ruostekoski et al. 



Q(X+iY) [Fig2b, Ruostekoski et al.] 




Q(X+iY) [Fig3b, Ruostekoski et al.] 




[Fig4a, Ruostekoski et al.] 
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